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ABSTRACT 

We consider fifty transiting short-period giant planets for which eclipse depths have been mea¬ 
sured at multiple infrared wavelengths. The aggregate dayside emission spectrum of these 
planets exhibits no molecular features, nor is brightness temperature greater in the near- 
infrared. We combine brightness temperatures at various infrared wavelengths to estimate 
the dayside effective temperature of each planet. We find that dayside temperatures are pro¬ 
portional to irradiation temperatures, indicating modest Bond albedo and no internal energy 
sources, plus weak evidence that dayside temperatures of the hottest planets are dispropor¬ 
tionately high. We place joint constraints on Bond albedo, Ab* and day-to-night transport 
efficiency, e, for six planets by combining thermal eclipse and phase variation measurements 
(HD 149026b, HD 189733b, HD 209458b, WASP-12b, WASP-18b, and WASP-43b). We con¬ 
firm that planets with high irradiation temperatures have low heat transport efficiency, and that 
WASP-43b has inexplicably poor transport; these results are statistically significant even if the 
precision of single-eclipse measurements has been overstated by a factor of three. Lastly, we 
attempt to break the Ab~£ degeneracy for nine planets with both thermal and optical eclipse 
observations, but no thermal phase measurements. We find a systematic offset between Bond 
albedos inferred from thermal phase variations (Ab ~ 0.35) and geometric albedos extracted 
from visible light measurements ( A g « 0.1). These observations can be reconciled if most 
hot Jupiters have clouds that reflect 30-50 per cent of incident near-infrared radiation, as well 
as optical absorbers in the cloud particles or above the cloud deck. 

Key words: methods: statistical - planets and satellites: atmospheres - infrared: planetary 
systems. 


1 INTRODUCTION 

Mature planets on short-period orbits have energy budgets domi¬ 
nated by incoming radiation rather than internal heat. Their atmo¬ 
spheric temperatures are therefore a function of both the absorption 
of incident stellar energy and its transport before re-emission into 
space. 

Absorbed energy is solely a matter of incident stellar flux and 
the planet's Bond albedo, Ab . Ironically, it is difficult to constrain 
Bond albedo (the fraction of stellar energy that is reflected) through 
observations of reflected light. In order to convert an optical geo¬ 
metric albedo (visible light reflected towards the illuminating star) 
into a Bond albedo, one must make assumptions about a planet’s re¬ 
flectance spectrum, scattering phase function, and spatial inhomo¬ 
geneity ( |Hanel|2003[ l. Bond albedo is more readily obtained from 
thermal measurements via energy balance. 

Heat transport is more complicated, but tends to move energy 
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from hot to cold: vertically upwards, from equator to pole (for plan¬ 
ets with small obliquity), from summer hemisphere to winter hemi¬ 
sphere (for planets with non-zero obliquity), and from day to night 
(for planets with slow rotation). Due to strong tides, short-period 
planets are expected to have zero obliquity and slow rotation, and 
most are on circular orbits. Hot Jupiters on circular orbits are fur¬ 
ther expected to be tidally-locked, with one side permanently fac¬ 
ing the host star and the other forever dark. As such, the atmosphere 
tends to transport heat from the dayside to the nightside, and from 
equatorial regions to the poles (for a recent review of hot Jupiter 
atmospheric dynamics, see jHeng & Showman|2014) . 

The atmospheric temperature of a planet is generally a func¬ 
tion of four variables: longitude, latitude, pressure (or height), 
and time. The time-dependence can usually be neglected for hot 
Jupiters because they exhibit minimal weather jAgol et al. |2010[ 
| Knutson et al.|201 l[[Wong et al.|2014) . Moreover, hot Jupiters on 
circular orbits are expected to have a 3D fixed temperature struc¬ 
ture with respect to the sub-stellar location, regardless of whether 
they are tidally locked {Rauscher & Kempton|2014j , Showman et al. 
|2014j>. This motivates using a star-based coordinate system with the 
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prime meridian facing the star, and allows us to use orbital phase 
as a proxy for longitude i jC'owan & Agol||2008fr . The latitudinal 
temperature-dependence of a hot Jupiter is inaccessible unless one 
can measure higher-order phase modulation i Cowan et al.|2013| or 
utilize occultation mapping i jMajeau et al.|2(F Ope Wit et aL[2012j l. 
Finally, the vertical temperature structure is in principle accessible 
via emission spectroscopy: wavelengths at which the atmosphere is 
relatively transparent will probe deeper layers, and vice versa. 

Multi-wavelength thermal phase variations of a hot Jupiter on 
a circular orbit therefore amount to brightness temperature mea¬ 
surements as a function of orbital phase and wavelength (e.g., 
Stevenson et a l.|2014[ l. If one is solely interested in the global prop¬ 
erties of the planet—namely Bond albedo and day-to-night heat 
transport—then one can further simplify the problem by combining 
brightness temperatures at each orbital phase to obtain a bolomet- 
ric flux, and hence an effective temperature at that phase. Note that 
Solar System planets tend to have effective temperatures that are 
either uniform from any vantage point, or which vary based on the 
latitude of the observer due to imperfect poleward heat transport. 
For short-period exoplanets, on the other hand, the principal tem¬ 
perature gradient is between day and night, and dayside effective 
temperatures are often hundreds to thousands of degrees greater 
than their nightside counterparts. The final simplification we make 
is therefore to treat the planet as two horizontally isothermal hemi¬ 
spheres: a dayside and a nightside. The effective temperatures of 
each hemisphere are simply the weighted mean of the atmospheric 
temperatures on that side of the planet. 


Table 1 . Short-period giant planets with a minimum of two published 
eclipse observations at infrared wavelengths (non-detections are not in¬ 
cluded.) 


Infrared Multi-Eclipse Planets 


CoRoT- lb 
CoRoT-2b 


GJ 436b 
HAT-P-lb 
HAT-P-2b 
HAT-P-3b 
HAT-P-4b 
HAT-P-6b 
HAT-P-7b 
HAT-P-8b 
HAT-P-12b 
HAT-P-19b 
HAT-P-20b 
HAT-P-23b 
HAT-P-3 2b 
HD 149026b 
HD 189733b 
HD 209458b 
KELT-lb 
Kepler-5b 
Kepler-6b 
Kepler- 12b 
Kepler-13 Ab 
Kepler- 17b 
TrES-lb 


TrES-2b 
TrES-3b 
TrES-4b 
WASP-lb 
WASP-2b 
WASP-3b 
WASP-4b 
WASP-5b 
WASP-6b 
WASP-8b 
WASP-12b 
WASP-14b 
WASP-17b 
WASP-18b 
WASP-19b 
WASP-24b 
WASP-26b 
WASP-33b 
WASP-39b 
WASP-43b 
WASP-48b 
XO-lb 
XO-2b 
XO-3b 
XO-4b 


1.1 Previous Work 

|Cowan & A gol (201 lb) used broadband infrared eclipse measure¬ 
ments of 24 hot Jupiters to demonstrate that they have generally 
low Bond albedos (Ab < 0.5), and that the hottest planets have 
extremely low albedos and/or poor day-night heat transport effi¬ 
ciency, £. 

It is possible to break the albedo-transport degeneracy by com¬ 
bining dayside thermal constraints with measurements of either 
nightside thermal emission or dayside reflected light. |Cowan et al.] 
< |2007} ) used an 8.0 pm eclipse depth and phase amplitude from the 
Spitzer Space Telescope, combined with an optical eclipse mea¬ 
surement from the MOST satellite, to constrain the energy budget 
of HD 209458b; they inferred a small Bond albedo (absorption of 
almost all light that shines on it) and a high day-night heat transport 
efficiency (nightside not much cooler than the dayside). 

|Cowan et al. 1(2012} used 8 infrared dayside and 2 mid-infrared 
nightside measurements to constrain the albedo and recirculation 
of WASP-12b; they found the planet has a modest Bond albedo 
(~ 0.25) and low heat transport efficiency (< 0.1). 

Stevenson et al.| ( |2014[ l used phase-resolved emission spec¬ 
troscopy taken with WFC3 from the Hubble Space Telescope to 
map the atmospheric thermal structure of WASP-43b, finding low 
Bond albedo (0.06-0.25) and no heat redistribution (for recent re¬ 
views of exoplanet atmospheric observations, please see Burrows 
|20l4| |B aileyj20T4| > . 

Our work is organized as follows: in Section [2] we use pub¬ 
lished eclipse depths at multiple infrared wavelengths to infer ef¬ 
fective dayside temperatures for fifty planets, with more than twice 
the data as Cowan & Agol (2011b). In Section [3] we consider the 
subset of planets for which we can break the albedo-recirculation 
degeneracy. We first tackle the six planets with thermal measure¬ 
ments of both eclipses and phase variations (Section[TTJ, then the 
nine planets for which reflected light measurements are available 


in addition to dayside thermal constraints (Section[T2]l. We discuss 
our results in Section [4] 


2 DAYSIDE ENERGY BUDGET 


Inferring the effective dayside temperature of a planet requires 
combining eclipse depths at thermal wavelengths, which we de¬ 
fine as those longward of 0.8 pm. We have updated the sample 
from Hansen et al.'|l j2014| l, now including fifty planets with a min¬ 
imum of two published infrared eclipse measurements (additions 
include HAT-P-19b, HAT-P-20b, HAT-P-32b, WASP-6b, WASP- 
26b, WASP-39b; |Mahtani et aTpOO] |Deming et aL|20l4l [Zhao| 
|et al.|2014| Kammer et al. in prep.). Our data predominantly consist 
of broadband photometry, but we also include spectroscopic emis¬ 
sion measurements when they are at complementary wavelengths 


(e.g. Ranjan et al.|2014||Wilkins et al.|2014||Crouzet et al.|2014[ 

IStevenson et al. 12014b. The planets from our sample are listed in 

Table □ 


2.1 Brightness Temperatures 

Thermal emission at different wavelengths originates from differ¬ 
ent layers in the planet’s atmosphere, which have different temper¬ 
atures. One can define a brightness temperature at each observed 
wavelength, Ti,(A): this is the temperature that a blackbody must 
have in order to emit at the same intensity as the planet. 

Inverting the Planck function, we obtain the following expres¬ 
sion for brightness temperature ( jCowan & Agol|201 lb} : 

he/XkT* - v -i - 1 

log ( 1 + 


U A) = 


he 


A k 


- 1 


■0(A) 


( 1 ) 


where t/>(A) is the relative intensity of the planet to that of its host 
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Combined Temperature Spectrum 




Figure 1. Average broadband emission spectrum for fifty short-period gi¬ 
ant planets (blue), plotted with the emission spectra of individual planets 
(gray). The spectrum of each planet is normalized to its mean brightness 
temperature. The aggregate spectrum is the median normalized brightness 
temperature at each wavelength, where the uncertainty bars denote the un¬ 
certainty on the mean (as opposed to the standard variation of the spectra at 
that wavelength.) 


star and T * is the stellar effective temperature, meaning we treat the 
star as a blackbody. For dayside measurements, ^>(A) is the ratio of 
eclipse depth to transit depth, <5 ec i/5tr, while for nightside measure¬ 
ments it is the ratio of nigthside flux to transit depth, (5 e ci — 5 V ar)/<5tr, 
where <5 var is the phase variation amplitude. Published data there¬ 
fore allow us to compute dayside (and, when appropriate, night- 
side) brightness temperatures for each waveband in which a planet 
has been observed. 


2.2 Aggregate Emission Spectrum 

The broadband emission spectra of most individual planets are con¬ 
sistent with isothermal atmospheres (Hansen et al. |2014) . It is pos¬ 
sible, however, to construct an aggregate emission spectrum of all 
fifty planets in the hopes of revealing molecular absorption features 
too faint to detect in any individual planet’s spectrum. 

If some planets have temperature inversions while others do 
not, this type of averaging could actually wash out molecular sig¬ 
natures, which would appear in absorption for some planets and 
emission for others. However, the first and most statistically sig¬ 
nificant case of a hot Jupiter temperature inversion i jKnutson et al.| 
|2008| has not been borne out by new measurements nor reanal¬ 
ysis of the originals (Zellem et al.|[2QI4[ [Diamond- Lowe et al.| 
|2014|[Schwarz et al.|20 1 5} . Moreover, a systematic study of Spitzer 
eclipse measurements found that they have not been as accurate as 
advertised ( [Hansen et al.|2014] ), suggesting the temperature inver¬ 
sions reported in most hot Jupiter atmospheres may simply be due 
to confirmation bias. 

The aggregate emission spectrum for the fifty hot Jupiters is 
shown in Figure |T] We normalize the brightness temperature spec¬ 
trum of each planet in our sample, then determine the median and 
uncertainty on the mean at each wavelength. This “stacking” is only 
useful, however, at wavelengths for which there are observations of 
many planets (currently 1.15, 1.65, 2.25, 3.6, 4.5, 5.8, and 8.0 mi¬ 
crons). 

There are no significant features in the average spectrum, not 
even the trend towards higher brightness temperatures at shorter 


Figure 2. Dayside effective temperature versus irradiation temperature for 
all giant planets with multiple dayside infrared eclipses, estimated by Monte 
Carlo using hybrid EWM-PWM calculation and inflating observational un¬ 
certainties by / s y S = 3 where applicable {Hansen et al.|20f4} . Uncertainty 
bars for both temperatures are shown, while dot size is proportional to the 
fraction of expected planetary emission that falls within observed wave¬ 
bands. Red symbols denote eccentric planets (e > 0.1). The solid, dashed, 
and dotted lines correspond respectively to maximum dayside tempera¬ 
ture, Td = (2/3) 1/,4 To, uniform dayside with zero nightside temperature, 
Td = (1/2) 1/,4 To, and equilibrium temperature, Td = (1 /4) 1 ’ ’’To. The 
trend line is shown in green; it suggests that hotter planets have dispropor¬ 
tionately hot daysides. 


wavelengths reported by Burrows (20T4). It is worth noting that 
Figure 4 of Burrows (2 014} used data from fewer planets, and was 
normalized differently: the equilibrium temperature of each planet 
was divided out, rather than its actual dayside effective temperature. 
Recall that a planet’s equilibrium temperature is what one would 
expect for a planet with zero Bond albedo and uniform tempera¬ 
ture; it is merely a convenient theoretical quantity proportional to 
the irradiation temperature, To, that we utilize in this work. The 
dayside effective temperature, on the other hand, is the weighted 
mean brightness temperature of the planet’s dayside, as described 
in Section [273] Most hot Jupiters have dayside effective tempera¬ 
tures greater than their equilibrium temperature (dotted line in Fig¬ 
ure [2} due to imperfect day-night heat transport. As a result, the 
dayside of a hot Jupiter emits somewhat more in the mid-IR—and 
considerably more in the near-IR—than one would predict based 
on its equilibrium temperature. In any case, we agree with Burrows 
( |2014) > that there are no signs of molecular absorption features in 
the aggregate spectrum. Since molecules are undoubtedly present 
in the atmospheres of exoplanets, we conclude that their absorp¬ 
tion features are being muted by vertically isothermal atmospheres, 
optically thick cloud, or both. 


2.3 Effective Temperatures 

While brightness temperatures of brown dwarfs are strongly 
wavelength-dependent in the near-infrared ( [Faherty et al.|[20l4| 
|Biller et al.||2013} , the external, asymmetric heating experienced 
by hot Jupiters produces dayside atmospheres that are relatively 
isothermal in the vertical direction (Fortney et ah| |2006| >. This 
results in relatively featureless dayside emission spectra, which 
are amenable to model-independent estimates of bolometric flux 
and hence effective dayside temperature. If the nightsides of hot 
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Jupiters have greater vertical temperature structure, then nightside 
effective temperature estimates are less reliable. 

There is no universal way to derive the effective temperature 
of a planet from a collection of brightness temperatures. We there¬ 
fore consider two methods with different physical motivations. In 
the first method, each brightness temperature is weighted by the in¬ 
verse square of its respective uncertainty (w, = 1/af), so eclipse 
depths with small relative uncertainties contribute more to the in¬ 
ferred effective temperature. We call this the error-weighted mean 
(EWM) effective temperature. This method has the advantage of 
being robust to occasional outlier eclipse depths, but implicitly as¬ 
sumes that short-period planets have Planck-like broadband spec¬ 
tra. 

The second method, which we only apply to dayside measure¬ 
ments, weighs the brightness temperatures by the expected inte¬ 
grated power in that waveband: uii = Pi = f* 2 3 B(T est , \)d\, 
where _B(T es t, A) is the Planck function for the planet's estimated 
dayside effective temperature. To bypass an iterative solution, we 
adopt Test = (l) 1 ^ 4 ^), where To is the irradiation temperature: 
To = T,^R/ / a where is the stellar radius and a is the semi¬ 
major axis. This is an excellent match (the dashed line in Figure|2j 
to the actual dayside effective temperature of most planets in our 
sample. We call this the power-weighted mean (PWM); it is iden¬ 
tical in spirit to the linear interpolation method of |Cowan & Agolj 
l |2011b) , but is easier to implement and runs faster. The PWM gives 

more weight to measurements near the peak of the planet’s Planck 
function and should, in the limit of high spectral coverage, pro¬ 
duce accurate effective temperatures even if planets have broad¬ 
band spectral features. 

Note that both EWM and PWM are biased in favor of broad¬ 
band measurements: these observations tend to have smaller uncer¬ 
tainties, and they capture more of the planet’s expected blackbody 
emission. The two methods produce generally consistent effective 
temperature estimates, which is a testament to the fact that most 
current dayside emission spectra are approximately Planck-like. 

We use a 10 4 -step Monte Carlo analysis to estimate uncertain¬ 
ties in dayside effective temperatures. At each step in the Monte 
Carlo, we randomly vary the stellar effective temperature, transit 
depth, eclipse depth, and scaled semi-major axis, a* = a/R f , 
according to their uncertainties. We use the published uncertain¬ 
ties for all of the above, except for single-epoch broadband eclipse 
depths where we inflate the published uncertainty by the factor 
/ sys = 3 (such measurements have historically been less accurate 
than advertised; |Hansen et al.|20l4) . We also compute each planet’s 
irradiation and brightness temperatures (following Equation^. We 
then estimate each planet’s dayside effective temperature; to hedge 
our bets, we use the EWM for half of the MC steps, and the PWM 
for the other half. The resulting relationship between To and T,/ is 
shown in Figure [2] The median property is plotted, and uncertainty 
bars denote standard deviation from the MC. 


3 GLOBAL ENERGY BUDGET 

Based on dayside effective temperatures alone, one cannot simulta¬ 
neously specify Bond albedo and heat recirculation efficiency. This 
degeneracy can be broken by supplementing thermal eclipses with 
one of two measurement types: phase variations at infrared wave¬ 
lengths, or eclipse depths at visible wavelengths. Table [2] lists the 
published data for the fifteen planets which fall into one or both 
of these categories: observed wavelength and bandwidth, eclipse 
depths, and phase amplitudes. Cyan-colored entries are measure- 



Dayside Temperature Ratio 


Reflection Ratio 



0 0.2 0.4 0.6 0.8 1 

Geometric Albedo 



Figure 3. Dayside temperature (top left), reflection (top right), and dayside 
flux for CoRoT-2b at 1.4 fim (bottom), shown as a function of albedo and 
recirculation (assuming the geometric and Bond albedos are equal to one 
another; see Section [43] for caveats.) At high albedo, the NIR dayside flux 
roughly parallels reflected starlight, while at low albedo, the dayside flux is 
mostly thermal emission and hence depends on day-night heat transport. 


ments exempt from the / sys = 3 uncertainty inflation ofQT ansen 
|etal1 ( |20T4) . 


3.1 Full-Orbit Thermal Measurement 


The first way to resolve the degeneracy between Bond albedo and 
heat recirculation is by combining thermal eclipse and phase mea¬ 
surements to infer the planet’s nightside effective temperature. This 
requires phase variations at thermal wavelengths, again defined 
as longward of 0.8 ;im. Such phase observations are more time- 
intensive than eclipses, and therefore less widely available. There 
are six planets with at least one published phase measurement: 


HD 149026b, HD 189733b, HD 209458b, WASP-12b, WASP-18b, 
and WASP-43b. As in Section [2] we include band-integrated spec¬ 
troscopy when it complements photometric observations (only one 
case at the moment: Stevenson et al.|2014) . For non-detections, an 
ncr upper limit of a is assumed to have a value and uncertainty 
of a/2 and a/(2n) respectively. Table [2] shows the data for this 
sample. 


Observational references are as follows: HD 149026b: Steven- 

son et al. 

(2012); Knutson et al. (2009); HD 189733b: Evans 

et al.; (20 

3); Crouzet et al. (2014); Barnes et al. (2007); Knut- 

son et al. i 

2012); Charbonneau et al. (2008); Todorov et al. (2014); 

Agol et a 

. (2010); Deming et a 

. (2006); Knutson et al. 

(2009 k 

HD 209458b: Rowe et al.[f2008); 

Richardson et al. (2003); 

Cnutson 

et al. (20081; Zellem et al. (2014); 

Cowan et al. (2007); WASP-12b: 

Fohring et al. (2013); Cro 1 et al. 

(2014); Crossfield et al. 

(20121; 

Stevenson et al. (2014); Swain et al. (2013;); Cowan et al. 

(20121; 


43b:|Stevenson et al.|j20 1 4); Wang et al. 1 2013); Zhou et al.|(|2014); 
IBlecic et al,| i |2014) . 
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Balancing the Energy Budget 5 


Table 2. Eclipse depths and phase amplitudes for our restricted planetary samples where degeneracy between albedo and heat recirculation can be 
resolved. Observations are denoted by central wavelength and bandwidth; measurements in cyan are exempt from the uncertainty inflation / sys = 3 
of |Hansen et al.|{2014) . 


Planet Wavelength (jzm) Eclipse Depth Phase Amplitude - Planet Wavelength (/xm) Eclipse Depth Phase Amplitude 


C 0 R 0 T-lb 

0.60(42) 

1.6(6) 


10- 4 



Kepler-7b 

0.65(40) 

3.870(835) 


10 -5 

4.8(13) 

X 10 -5 


1.65(25) 

1.45(49) 

X 

10 -3 




3.60(75) 

3.08(103) 

X 

10- 4 




2.10(2) 

2.8(5) 

X 

10 -3 




4.5(10) 

5.05(168) 

X 

10- 4 




2.15(32) 

3.190(405) 

X 

10 -3 



Kepler-13 Ab 

0.65(40) 

1.720(18) 

X 

10 -4 

1.5200(105) 

X 10 -4 


3.60(75) 

4.15(42) 

X 

10 -3 




2.15(32) 

1.22(51) 

X 

10 -3 




4.5(10) 

4.82(42) 

X 

10 -3 




3.60(75) 

1.56(31) 

X 

10 -3 



CoRoT-2b 

0.60(42) 

6(2) 

X 

10 -5 




4.5(10) 

2.22(23) 

X 

10 -3 




1.4(6) 

3.95(57) 

X 

10 -4 



TrES-2b 

0.65(40) 

7.7(18) 


10“ 6 

4.10(105) 

X 10 -6 


1.65(25) 

8.50(283) 

X 

10- 4 




2.15(32) 

6.2(12) 

X 

10- 4 




2.15(32) 

1.6(9) 

X 

10 -3 




3.60(75) 

1.27(21) 

X 

10 -3 




3.60(75) 

3.55(20) 

X 

10 -3 




4.5(10) 

2.30(24) 

X 

10 -3 




4.5(10) 

4.75(19) 

X 

10 -3 




5.8(14) 

1.99(54) 

X 

10 -3 




8.0(29) 

4.09(80) 

X 

10 -3 




8.0(29) 

3.59(60) 

X 

10 -3 



HAT-P-7b 

0.65(40) 

7.12(15) 

X 

10 -5 

7.33(27) 

X 10 -5 

WASP-12b 

0.90(15) 

1.360(136) 

X 

10 -3 




3.60(75) 

9.8(17) 

X 

10- 4 




1.04(12) 

1.09(14) 

X 

10 -3 




4.5(10) 

1.59(22) 

X 

10 -3 




1.25(16) 

1.39(30) 

X 

10 -3 




5.8(14) 

2.45(31) 

X 

10 -3 




1.38(55) 

1.580(39) 

X 

10 -3 




8.0(29) 

2.25(52) 

X 

10 -3 




1.4(6) 

1.740(17) 

X 

10 -3 



HD 149026b 

3.60(75) 

4.0(3) 

X 

10 -4 




1.65(25) 

1.91(20) 

X 

10 -3 




4.5(10) 

3.4(6) 

X 

10 -4 

2.23(58) 

X 10 -4 


2.15(32) 

2.96(14) 

X 

10 -3 




5.8(14) 

4.4(10) 

X 

10- 4 




2.220(34) 

3.01(46) 

X 

10 -3 




8.0(29) 

3.7(8) 

X 

10 -4 

2.3(7) 

X 10 -4 


2.320(27) 

4.5(6) 

X 

10 -3 




16(6) 

8.5(32) 

X 

10 -4 




3.60(75) 

4.19(44) 


10 -3 

3.20(33) 

X 10 -3 

HD 189733b 

0.37(16) 

1.26(37) 

X 

10 -4 




4.5(10) 

4.29(33) 

X 

10 -3 

3.92(16) 

X 10 -3 


0.51(12) 

1(34) 

X 

10 -6 




5.8(14) 

6.96(60) 

X 

10 -3 




1.4(6) 

9.6(39) 

X 

10 -5 




8.0(29) 

6.96(96) 

X 

10 -3 




2.15(32) 

2(2) 

X 

10- 4 



WASP-18b 

3.60(75) 

3.04(26) 

X 

10 -3 

2.96(11) 

X 10 -3 


3.60(75) 

1.47(4) 


10 -3 

1.240(61) 

X 10 -3 


4.5(10) 

3.79(21) 


10 —3 

3.66(9) 

X 10 -3 


4.5(10) 

1.790(38) 

X 

10 -3 

9.82(89) 

X 10 -4 


5.8(14) 

3.7(3) 

X 

10 -3 




5.8(14) 

3.10(34) 

X 

10 -3 




8.0(29) 

4.1(2) 

X 

10 -3 




6.45(210) 

2.200(62) 

X 

10 -3 



WASP-19b 

0.685(530) 

3.9(19) 

X 

10~ 4 




8.0(29) 

3.44(36) 

X 

10 -3 




0.91(20) 

8.0(29) 

X 

10- 4 




10.5(60) 

3.560(67) 


10 -3 




1.29(8) 

8.3(39) 

X 

10 -4 




16(5) 

5.51(30) 

X 

10 -3 




1.6(4) 

1.86(14) 

X 

10 -3 




24(9) 

5.36(27) 

X 

10 -3 

1.3(3) 

X 10 -3 


1.65(25) 

2.76(44) 

X 

10 -3 



HD 209458b 

0.5(3) 

7(9) 

X 

10 -6 




2.15(32) 

2.87(20) 

X 

10 -3 




2.15(32) 

1.5(15) 

X 

10- 4 




3.60(75) 

4.83(25) 

X 

10 -3 




3.60(75) 

9.4(9) 

X 

10 -4 




4.5(10) 

5.72(30) 

X 

10 -3 




4.5(10) 

1.3900(705; 


10 -3 

1.090(115) 

X 10 -3 


5.8(14) 

6.5(11) 

X 

10 -3 




5.8(14) 

3.01(43) 

X 

10 -3 




8.0(29) 

7.3(12) 

X 

10 -3 




8.0(29) 

2.40(26) 

X 

10 -3 

7.50(375) 

X 10 -4 

WASP-43b 

1.4(6) 

4.61(5) 

X 

10 -4 

4.68(4) 

X 10 -4 


24(9) 

3.38(26) 

X 

10 -3 




1.65(25) 

1.03(17) 

X 

10 -3 



Kepler-5b 

0.65(40) 

1.86(36) 

X 

10 -5 

1.93(58) 

X 10 -5 


2.15(32) 

1.81(27) 

X 

10 -3 




3.60(75) 

1.03(17) 

X 

10 -3 




3.60(75) 

3.47(13) 

X 

10 -3 




4.5(10) 

1.07(15) 

X 

10 -3 




4.5(10) 

3.82(15) 

X 

10 -3 



Kepler-6b 

0.65(40) 

1.11(40) 

X 

10 -5 

1.72(43) 

X 10 -5 









3.60(75) 

6.9(27) 

X 

10- 4 











4.5(10) 

1.51(19) 

X 

10 -3 











3.1.1 Reflected Infrared Light 

The light emanating from a planet’s dayside is a combination of 
thermal emission and reflected starlight. We plot an example for 
CoRoT-2b at 1.4 /im in the bottom panel of Figure [3] assuming 
the geometric and Bond albedos are equal (though as described in 
Section [43] the proper conversion is more involved.) If the planet 
has low albedo in this scenario, then the 1.4 /mi flux is almost en¬ 
tirely thermal emission and depends primarily on day-night heat 
transport. In the high albedo limit, on the other hand, the NIR flux 
is primarily reflected light and so varies linearly with the geomet¬ 
ric albedo. In other words, even eclipse measurements at wave¬ 
lengths greater than 0.8 /jm are potentially contaminated by re¬ 
flected starlight. 

Furthermore, Figure [4] shows the reflected light contribution 
to dayside flux as a function of wavelength for a hypothetical gray 
planet (with temperature limits derived from Equation [2] see Sec¬ 
tion |3.1.2| ) Reflected light dominates at ultraviolet wavelengths 
as expected, but its prevalence continues well through the near- 
infrared. For reasonable system parameters, reflected light con¬ 
tributes > 10 per cent of the NIR flux (this reflected contribution 
goes up if T* is increased, or if Td or a* are decreased.) Even if 
molecular absorption depresses the reflectance in certain bands, it 
is likely that eclipse measurements in NIR water opacity windows 


(J, H, and K ) are contaminated by reflected light at the > 10 per 
cent level. 


3.1.2 Confidence Regions 

In estimating effective dayside and nightside temperatures for each 
planet, we use Monte Carlo simulations with 5000 steps to propa¬ 
gate uncertainties. We assume uncertainties on observed quantities 
to be Gaussian and symmetrical; when asymmetrical uncertainty 
bars are reported, we adopt their mean. The planetary irradiation 
temperature is first computed as described in Section [273] Bright¬ 
ness temperatures are calculated for each appropriate measurement 
front Table [2] using Equation [Tj propagating uncertainty on stellar 
effective temperature, eclipse depth, phase amplitude, and transit 
depth. For observations not listed in cyan (isolated eclipses, or par¬ 
tial phase curves), we conservatively inflate the published uncer¬ 
tainty by the factor f sys = 3 (Hansen et al.|2014| >. Reflected light 
contributions are subtracted from all dayside eclipse depths—using 
planetary radius, R p , and semi-major axis—assuming infrared geo¬ 
metric albedos to be normally distributed with mean 0.07 and width 
0.01 (this is the distribution of uncorrected optical geometric albedo 
values described in Section [4~5] i. In cases where a brightness tem¬ 
perature is calculated as 0 K for all MC steps, we assume 100 K 
uncertainty in subsequent propagations. We then compute the ef- 
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Wavelength (ftm) 


A g 

■ 0.02 
■ 0.1 
■ 0.25 
□ 0.5 


Figure 4. Contribution of planetary reflected light as a function of wave¬ 
length for different geometric albedos (lig hter = higher), assuming black- 
body emission and q = | (see Section 4. 5|. This example planetary system 
assumes T* = 6100 K and a* = 4.8, both weighted means of the fif¬ 
teen planets in Tabled] Albedo regions are bounded by dayside temperature 
limits using Equation^] solid lines denote no recirculation (e = 0), dashed 
lines denote perfect recirculation (e = 1). The vertical dashed line indicates 
our cutoff of 0.8 fim between reflected light and thermal emission. 


fective dayside and nightside temperatures using the EWM of the 
corresponding brightness temperatures (as this requires no a priori 
temperature assumption and produces similar values to the PWM.) 

Our parameterization of recirculation neglects any treatment 
of poleward heat transport. The most extreme meridional temper¬ 
ature profiles are either uniform in the North-South direction (per¬ 
fect poleward transport) or T oc cosi 9, where 9 is latitude (no 
poleward transport). The difference in effective temperature seen 
by an equatorial observer is (l/4) 1,/4 Tb versus (8/37T 2 ) 14 To, a 
1 per cent discrepancy. We incorporate this worst-case systematic 
uncertainty in quadrature for all effective temperature estimates. 

Understandably, the number of brightness temperature mea¬ 
surements at distinct wavelengths for a planet affects the accuracy 
of the effective temperature estimate. In a Monte Carlo analysis 
using J.J. Fortney atmospheric models, |Cowan & Agol|(20 1 lb[ l es¬ 
timated systematic errors of 7.6 per cent in effective temperature 
when only a single observation was used (note that we only con¬ 
sider planets with at least two measurements), down to approxi¬ 
mately 2.5 per cent for four or more measurements. We conser¬ 
vatively adopt a similar sliding scale of 8 per cent down to 3 per 
cent systematic uncertainty in effective temperature over the same 
observation number range, again added in quadrature. 

Once we have dayside and nightside effective temperatures— 
and realistic uncertainties—for the six exoplanets, it is possible to 
infer each planet’s Bond albedo and day-night heat transport effi¬ 
ciency using the parameterization of Cowan & Agol ( 201 lb[ >: 

T^Toil-An) 1 '*^- ' , ( 2 ) 

and 

T n = T 0 (l~A B ) 1/i (^y / \ (3) 

where both Ab and e can take values between zero and unity. 

We create \ 2 surfaces for each planet based on our estimated 
dayside and nightside effective temperatures and using Equations[2] 
and H] We calculate y 2 on a 101 x 101 grid in A B and e, then in¬ 
terpolate the intermediate values. The la, 2a, and 3a confidence 
intervals are defined as A\ 2 = {1,4,9} respectively above the 
minimum, \7n in , where Xmin ~ 0 for most planets because we have 



Albedo 

Figure 5. Composite la confidence regions for thermal observation plan¬ 
ets, as calculated from the error-weighted mean dayside and nightside 
brightness temperatures. Here the horizontal axis measures Bond albedo. 
Bounding curve colors indicate irradiation temperature: red = warmer, pur¬ 
ple = cooler. The inflationary factor f sys = 3 is applied to infrared eclipse 
and phase uncertainties as noted in Table [2] but adopting the published 
eclipse uncertainties barely modifies the confidence intervals. 


two constraints and two model parameters. Generating the \ 2 sur¬ 
faces involves numerical integrations of Planck functions, which 
can be computationally intensive. We therefore create a database of 
relevant integrals; with 10 4 grid points tested per effective temper¬ 
ature, this database decreases computational time by more than 95 
per cent. 

We plot the la confidence intervals for the six exoplanets with 
full-orbit thermal measurements in Figure [5] Each planet is col¬ 
ored according to irradiation temperature, essentially the incident 
stellar flux. Since these planets have benefited from intensive ob¬ 
servational campaigns, omitting the / sys = 3 uncertainty inflation 
produces nearly identical confidence intervals. 


3.2 Geometric Albedo Measurement 


The alternative approach to resolving the albedo versus heat- 
transport degeneracy of thermal eclipses is to also consider eclipse 
measurements at optical wavelengths. For our purposes, optical 
eclipses are those shortward of 0.8 gm; these observations allow 
us to infer the planet’s optical geometric albedo. Our literature re¬ 
view finds nine planets with published eclipse depths at thermal 
and optical wavelengths, but lacking infrared phase variations (Ta¬ 
ble^. Note that HD 189733b and HD 209458b benefit from both 


infrared phases and visible eclipses. 


2009 Zhao et al. 2012; Gillon, M. et al. 2009; Rogers et al. 2013V 

Deming et al. 20111, CoRoT-2b (Alonso, R. et al. 2009 Wilkins 

et al. 2014 

|Alonso et al. 2010 [ Deming et al.|2011), HAT-P-7b (Es- 

teves et al. 

2014, Christiansen et al. 2010), Kepler-5b (Esteves et al. 

2014 

Desert et al. 2011), Kepler-6b (lEsteves et al. 2014, Desert 

et al. 

2011), Kepler-7b (Esteves et al. 2014, Demory et al. 2013), 


-2b 


( |Esteves et al. 2014HCroll et al.|2010| |ODonova net al.|2010| >, and 
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Wavelength (/rm) 


Td (K) 

— 500 

— 1000 

— 1500 

— 2000 

— 2500 

— 3000 


Figure 6. Flux ratio between blackbody and “Sum of Blackbodies” spectra 
for various effective dayside temperatures in 500 K increments. Curves 
are colored according to temperature: red = warmer, purple = cooler. The 
vertical dashed line indicates our chosen threshold wavelength, 0.8 ^zm, 
between reflected light and thermal emission. 


WASP-19b (Abe, L. et al.T20T3l|Zhou et al.|2013||Bean et al.|2013| 
|Zhou et al. 12014) . 


3.2.1 Thermal Contamination 

In order to extract a geometric albedo from an optical eclipse, we 
must correct the eclipse depth for thermal emission from the planet 
“leaking” into the visible band (Cowan & Ago! 12011b, Heng & 
|Demory|20 1 3} . In practice, one estimates the planet’s dayside ef¬ 
fective temperature and extrapolates this into the optical to account 
for thermal emission at visible wavelengths. However, this proce¬ 
dure is complicated by the fact that real hot Jupiters are vertically 
and horizontally inhomogeneous, so they emit at higher brightness 
temperatures in the optical than in the thermal infrared. 

Even if every location on a planet emits as a blackbody (BB), 
the resulting spectrum will not be a Planck curve. For a planet in the 
zero-albedo and zero-recirculation limit, the equilibrium tempera- 

i 

ture at any dayside location is described by T = To cos 4 7, where 
7 is the angle from the sub-stellar point (7 =3 5 at the termina¬ 
tor.) Each annulus of the dayside thus radiates at a different black¬ 
body temperature, and together they produce a “Sum of Blackbod¬ 
ies” (SoB) spectrum (this is analogous to the multicolor blackbody 
spectra used to model accretion disks; |Mitsuda et al.p984) . For 
fixed bolometric flux, BB and SoB spectra produce nearly identical 
flux at thermal wavelengths: the SoB is 1-2 per cent fainter than 
the BB. At optical wavelengths, however, a BB spectrum underes¬ 
timates the flux by a factor of a few, as seen in Figure[6] 

Moreover, the optical photosphere should be deeper and hotter 
than the infrared photosphere, in a cloud-free atmosphere ( |Allard| 
|et al.|200T]|Fortney et al.|2008||Cowan & Agol|201 la} . The combi¬ 
nation of these two effects is that a naive blackbody extrapolation 
from the infrared to the optical may underestimate thermal emis¬ 
sion by a factor of 3-10. In other words, while the hottest planets 
have the greatest thermal emission at optical wavelengths, some¬ 
what cooler planets have optical emission that is harder to estimate. 

Once an optical eclipse has been corrected for thermal con¬ 
tamination, the geometric albedo can be calculated using 



where is the reflected light eclipse depth. Note that geometric 
© 2015 RAS, MNRAS 00Q.[T|fT2l 



Albedo 


Figure 7. Composite 1 o confidence regions for eclipse-only planets as cal¬ 
culated from brightness temperatures (using EWM) and geometric albedos, 
shown with thermal observation planets for comparison. Here the horizontal 
axis measures different quantities: geometric albedo at visible wavelengths 
for eclipse-only planets (dashed lines). Bond albedo for thermal observation 
planets (solid lines). Bounding curve color follows Figure[5] The inflation¬ 
ary factor / s y s = 3 is applied to infrared eclipse and phase uncertainties as 
noted in Table[2] Adopting published uncertainties across the board results 
in similar confidence intervals. 


albedo is a function of wavelength, while hot Jupiter eclipse depths 
have typically only been measured in a single optical broadband. 


3.2.2 Confidence Regions 

Our dayside temperature analysis for planets in the eclipse-only 
sample is analogous to Section [3.1.2| and we again perform Monte 
Carlo simulation with 5000 steps for all uncertainty propagation. 
For optical eclipses, we optimize computation by first calculat¬ 
ing the uncertainty of each thermally-corrected geometric albedo, 
using Equation [4] and propagating uncertainties in EWM dayside 
effective temperature, stellar effective temperature, transit depth, 
eclipse depth, planetary radius, and semi-major axis. Our thermal 
correction uses equal contributions of BB and SoB spectra, assum¬ 
ing 10 per cent increase in dayside temperature to account for the 
vertical temperature profile effect noted in Section |3.2.1| To ac¬ 
knowledge variability with this effect, we add a 5 per cent system¬ 
atic uncertainty in quadrature to the calculated geometric albedo 
uncertainty. We then construct \ 2 surfaces for each planet as de¬ 
scribed in Section |3.1.2| For optical measurements, we recompute 
our thermal contamination correction at each \ 2 grid point to de¬ 
termine specific values of geometric albedo (since the underly¬ 
ing dayside temperature varies in the albedo-recirculation plane.) 
Note that HD 189733b has two distinct optical eclipses; we use 
the weighted mean of both corrected geometric albedos in our \ 2 
calculations for this planet. 

In Figure [ 7 ] we compare the la confidence intervals of the 
nine eclipse-only planets to those of the six planets with full-orbit 
thermal observations. Note that for optical eclipses we implicitly 
assume Ab = A g , but the actual conversion between geometric 
and Bond albedo is more complicated (Section |4.5fr . Regions are 
again colored by irradiation temperature, while sample group is 
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Table 3. Resulting parameters for all planets as calculated from brightness tempera¬ 
tures (via EWM) and geometric albedos, assuming applicable / sys = 3 uncertainty 
inflation. Thermal observation planets are listed first (with Bond albedos), followed 
by eclipse-only planets (with geometric albedos at visible wavelengths.) Low and 
high values are obtained from the confidence regions of Figures [5] and fit values 
are taken to be the grid location of Xmm- 


Planet 

Low 

Albedo 

Fit 

High 

Recirculation 

Low Fit High 

v 2 . 

Amin 

HD 149026b 

0.27 

0.41 

0.519 

0.045 

0.15 

0.353 

0.0005 

HD 189733b 

0.325 

0.37 

0.407 

0.536 

0.59 

0.648 

0.0051 

HD 209458b 

0.323 

0.43 

0.521 

0.283 

0.44 

0.606 

0.0013 

WASP-12b 

0.273 

0.37 

0.456 

0.024 

0.07 

0.172 

0.0028 

WASP-18b 

0 

0 

0.054 

0.007 

0.01 

0.032 

2.2571 

WASP-43b 

0.104 

0.29 

0.45 

0 

0 

0.024 

0.7581 

CoRoT-lb 

0 

0.01 

0.117 

0 

0 

0.232 

0.0002 

CoRoT-2b 

0.008 

0.07 

0.13 

0.246 

0.4 

0.545 

0.0005 

HAT-P-7b 

0 

0.04 

0.14 

0.472 

0.61 

0.745 

0.0012 

Kepler-5b 

0 

0.04 

0.107 

0.476 

0.67 

0.844 

0.001 

Kepler-6b 

0 

0.02 

0.082 

0.267 

0.51 

0.733 

0.005 

Kepler-7b 

0.258 

0.34 

0.426 

0.902 

1 

1 

0.9611 

Kepler-13 Ab 

0 

0.18 

0.341 

0.4 

0.58 

0.761 

0.0004 

TrES-2b 

0 

0.01 

0.056 

0.532 

0.67 

0.795 

0.0097 

WASP-19b 

0 

0.08 

0.191 

0.404 

0.56 

0.688 

0.0009 


denoted by the line style of bounding curve. Though the / sys un¬ 
certainty inflation is included for isolated thermal eclipses, taking 
the published uncertainties at face value produces nearly identical 
confidence intervals. 

Based solely on dayside effective temperatures (Figure[2]>, one 
might conclude that all planets have roughly the same Bond albedo 
and heat transport efficiency. Figure [7] dispels this notion at high 
significance. The lcr intervals for all the planets in Figure [7] are 
listed in Table [3] thermal observation planets first, eclipse-only 
planets second. All best-fit parameters are defined as the location 
°f A'min on the computed grids. 


4 DISCUSSION 

4.1 Sources of Error and Uncertainty 

It is worth summarizing the various sources of uncertainty and er¬ 
ror that we account for in order to produce Figure For thermal 
measurements, we first compute brightness temperatures, account¬ 
ing for the uncertainties on eclipse depth (inflated by a factor of 3 
if based on a single occultation; |Hansen et al.|2014) , transit depth, 
and stellar effective temperature, and also compute each planet’s 
irradiation temperature, accounting for uncertainty in stellar effec¬ 
tive temperature and scaled semi-major axis. At this stage, we also 
account for reflected light contamination (non-zero near-infrared 
geometric albedo). Nightside brightness temperatures are derived 
the same way, but additionally depend on the thermal phase am¬ 
plitude and its uncertainty. In converting brightness temperatures 
to effective temperatures, we account for unknown meridional heat 
transport and incomplete spectral coverage. The conversion from 
dayside brightness temperatures to dayside effective temperature 
is reasonable for hot Jupiters because of their relatively isother¬ 
mal vertical structure; the conversion may be more fraught for the 
nightsides of hot Jupiters. 

For optical eclipses, we first correct eclipse depths for thermal 
contamination, accounting for uncertainty on dayside temperature, 


transit depth, and stellar effective temperature. In this process we 
also account for vertical and horizontal temperature profiles that 
conspire to increase optical thermal emission. We next convert the 
reflected light eclipse depth to an optical geometric albedo, ac¬ 
counting for uncertainties in eclipse depth, planetary radius, and 
scaled semi-major axis. We assume Ag = A g for the purposes of 
constraining heat transport in Figure[7](see Section |4~5l for caveats), 
but this assumption in no way affects our inferred geometric albedo 
for these planets. 

Crucially, for every “correction” that we apply, we add ap¬ 
propriate uncertainty in our inference, either by randomly varying 
parameters in the Monte Carlo, or by adding systematic uncertainty 
in quadrature to the formal errors. Our inferences of heat transport, 
Bond albedo, and geometric albedo are therefore conservative. 


4.2 Dayside Temperatures 

The upward trend in dayside effective temperature with irradia¬ 
tion temperature in Figure [2] is unsurprising: one expects highly- 
irradiated planets to be hotter. The black lines in the plot can be 
thought of as limiting cases of either heat recirculation or Bond 
albedo. In the zero-albedo limit, the black lines correspond to e = 0 
(solid), e = 0.4 (dashed), and e = 1 (dotted). Alternatively, in 
the zero-recirculation limit, the black lines correspond to Ag =0 
(solid), Ag = 0.25 (dashed), and Ag = 0.625 (dotted). Therefore, 
planets that lie above the solid black line must have an internal en¬ 
ergy source, while planets lying below the dotted black line must 
have non-zero Bond albedo. 

We can also consider the qualitative claim from |Cowan &| 
|Agol| < |201 lbfr that Td increases disproportionately with To. We ig¬ 
nore planets with significantly eccentric orbits, as this complicates 
their energy budget (these planets are denoted in red in Figure [2]) 
With double the planets and more data per planet, we find that 
Td = —90(80) + 0.87(5)To-The x 2 P er datum of the fitted bend 
is 1.4 ± 0.4, which is a reasonable fit. This is consistent with—but 
does not strengthen—the claim that planets receiving more stellar 
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Table 4. Geometric albedo values, optical waveband starlight fraction, and Bond albedo limits (using Equations[6j{8| for planets with visual 
eclipses. Observations are denoted by central wavelength and bandwidth. Geometric albedo values are considered as follows: “uncorrected” 
uses Equation[4]assuming no thermal contamination; “simple correction” accounts for thermal emission with a BB spectrum at our error- 
weighted mean effective dayside temperature; “full correction” uses the method described in Section [3.2.2| (equal BB + SoB contributions 
with vertical temperature profile effect), again using our EWM dayside temperature. Uncertainties on starlight fraction are negligible and 
thus omitted. All Bond albedo calculations use the “full correction” values and assume the phase integral factor q = |. 


Planet 

Wavelength (^m) 

Uncorrected 

Geometric Albedo 
Simple Correction 

Full Correction 

fT 

Bond Albedo 

Min Gray High 

CoRoT- lb 

0.60(42) 

0.213(87) 

0.129(82) 

0.043(77) 

0.485 

0.026 

0.053 

0.284 

CoRoT-2b 

0.60(42) 

0.101(35) 

0.090(34) 

0.069(60) 

0.472 

0.04 

0.086 

0.305 

HAT-P-7b 

0.65(40) 

0.262(60) 

0.156(54) 

0.051(73) 

0.438 

0.028 

0.064 

0.309 

HD 189733b 

0.37(16) 

0.374(113) 

0.374(113) 

0.374(124) 

0.091 

0.043 

0.468 

0.497 


0.51(12) 

0.042(60) 

0.043(61) 

0.043(79) 

0.123 

0.007 

0.054 

0.445 

HD 209458b 

0.5(3) 

0.039(37) 

0.039(36) 

0.039(62) 

0.384 

0.019 

0.049 

0.327 

Kepler-5b 

0.65(40) 

0.107(22) 

0.079(25) 

0.040(58) 

0.439 

0.022 

0.05 

0.302 

Kepler-6b 

0.65(40) 

0.060(22) 

0.047(23) 

0.025(55) 

0.436 

0.014 

0.031 

0.296 

Kepler-7b 

0.65(40) 

0.314(73) 

0.315(72) 

0.313(88) 

0.44 

0.172 

0.392 

0.452 

Kepler-13 Ab 

0.65(40) 

0.468(31) 

0.318(67) 

0.175(113) 

0.405 

0.088 

0.218 

0.386 

TrES-2b 

0.65(40) 

0.030(7) 

0.022(8) 

0.007(51) 

0.439 

0.004 

0.008 

0.284 

WASP-19b 

0.685(530) 

0.248(120) 

0.177(113) 

0.099(108) 

0.537 

0.067 

0.124 

0.298 


flux generally have lower Bond albedo and/or less efficient heat 
transport. 


4.3 Thermal Phase Measurements 


Figure[5]shows a tendency towards lower recirculation efficiency as 
irradiation increases, in agreement with previous findings (Cowan 
& Agol||2011b| |Cowan et al.||2012| |Perez-Becker & Showman 
2013| >. The irradiation temperatures of these planets span approx¬ 
imately 2000 K, corresponding to e = 0.59 for HD 189733b 
at the cool end (T 0 = 1695 K) and e = 0.01 for WASP-18b 
(T 0 = 3387 K). The irradiation of WASP-12b is actually ~ 260 K 
higher than WASP-18b, but their recirculation probability distribu¬ 
tion functions overlap (Table[3]». 

The notable exception to the Tq-e trend is WASP-43b, with 
To = 1943 I\ but exhibiting virtually no heat transport (e = 0 
with Xmin — 0.758). Our recirculation value is in agreement with 
the redistribution factor of Stevenson et ah|( |2014| >, and our best- 
fit Bond albedo (Ab = 0.29) is also consistent at the la level. 
These parameters translate into a cold nightside temperature (nom¬ 
inally T n < 465 K), which we routinely find to be consistent with 
zero. Coupled hydrodynamic and radiative transfer simulations of 
this planet were able to reproduce its dayside—but not nightside— 
emission i jKataria et al.||2014} , so the poor heat transport of this 
planet remains a mystery. 

The thermal measurements for WASP-18b suggest 
Ab < 0.05 at 1 <t. The planet’s best-fit parameters would 

lie outside the plot to the left, which is indicative of either an 
internal heat source (identical to a negative Bond albedo in our 
parametrization) or underestimated observational uncertainties. 
Kepler-7b and WASP-43b also have Xmin well above zero, 
suggesting that either our model assumptions or the published 
uncertainties are incorrect. Each of these planets would benefit 
from more thermal eclipse and phase measurements in order 
to reduce the systematic uncertainty in dayside and nightside 
effective temperatures. 

HD 189733b and HD 209458b benefit from full-orbit ther¬ 
mal observations as well as optical eclipse measurements, allowing 
us to compare infrared-based Bond albedos to their optical geo¬ 
metric albedos. For HD 189733b we derive corrected geometric 


albedos of 0.37(12) at 0.37 fj ,m and 0.04(8) at 0.51 pm, in agree¬ 
ment with Evans et al.; (20T3j. This red-optical geometric albedo 
is much lower than our Bond albedo estimate of [0.33,0.41], For 
HD 209458b we obtain a corrected geometric albedo of 0.04(6) 
at 0.5 pm, which agrees with |Rowe et al.| (2008). However, our 
la interval for Bond albedo is [0.32,0.52], Therefore, the tentative 
conclusion based on these two planets is that their Bond albedos 
are considerably higher than their optical geometric albedos. 


4.4 Eclipse-Only Measurements 


Most of the eclipse-only planets have low geometric albedos: 
Ag < 0.2 (Figure [7j. As anticipated, Kepler-7b lies completely 
above this range ( [Demory et al.|20fT) , while the confidence region 
for Kepler-13 Ab extends to a geometric albedo of 0.34. The nine 
planets exhibit a wide variety of recirculation efficiencies, from 
CoRoT-lb (e w 0.1) to Kepler-7b (e w 0.95). Eclipse-only plan¬ 
ets with similar irradiation temperatures are found at different lo¬ 
cations on the e-axis, and we do not see evidence for a trend be¬ 
tween irradiation temperature and recirculation efficiency as with 
the thermal observation planets. This is not surprising, since day- 
side measurements offer minimal leverage for inferring the night¬ 
side temperature. Strictly speaking we only include these planets 
in Figure [7] by assuming that Ab = A g ; the actual comparison is 
more complex (Section [4.5fr . 

Geometric albedo analyses encompassing several planets from 
our sample have been previously conducted. We compare our re¬ 
sults in Table |4| to overlapping planets from |Heng & Demory| 
|20T3) : HAT-P-7b, Kepler-5b, Kepler-6b, Kepler-7b, and TrES- 
2b. Our “uncorrected” geometric albedos for all five planets show 
agreement within the stated confidence intervals. Esteves et al. 
< |20T4> and[Angerhausen et al. (2014) also consider these planets, 
in addition to Kepler-13Ab, under both zero and perfect heat redis¬ 
tribution. Our “full correction” geometric albedos for five of the six 
planets are in agreement with values from both studies obtained in 
the maximum equilibrium temperature hypothesis (i.e. hotter day- 
side temperatures implying greater thermal contamination of the 
optical eclipse). However, we find Kepler-13Ab to have dissimi¬ 
lar geometric albedo when using stellar parameters from |Shporer[ 
et al. 1 2014) with our greater thermal correction: 0.175(113) ver- 
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sus 0.404(55) and ~ 0, respectively. Note this includes stellar ef¬ 
fective temperature readjustment of Kepler-13A to 7650 ± 250 K, 
down from 8500 4= 400 K jSzabo et al.|20l l| i. 

4.5 Comparing Geometric to Bond Albedo 

In principle, an optical eclipse depth is related to a planet’s Bond 
albedo < |Rowe et al.||2006) . Indeed, for a range of Solar System 
planets and moons, the Bond albedo is roughly equal to the optical 
geometric albedo, albeit with a scatter of ±30 per cent. Given the 
possibility of inhomogeneous albedo, uncertainty in the scattering 
phase function, and unknown reflectance spectrum, it would be im¬ 
prudent to extrapolate this trend to hot Jupiters. Moreover, |Marley| 
|et al.|(l999[ l demonstrated that simply varying the incident stellar 
radiation can alter Bond albedo by a factor of four for identical 
planets. 

Our analysis of hot Jupiters suggests that their optical geomet¬ 
ric albedos are systematically lower than their Bond albedos. There 
are three possible explanations for this discrepancy: (1) we have 
over-corrected the thermal contamination at optical wavelengths, 
(2) we have systematically underestimated the effective tempera¬ 
tures for planets with full-orbit phase variations, or (3) the geomet¬ 
ric albedos of hot Jupiters are, in fact, systematically lower than 
their Bond albedos because of unexpected scattering phase func¬ 
tions and/or reflectance spectra. 

We address the first hypothesis by listing three different ge¬ 
ometric albedo calculations in Table [4] these differ in their treat¬ 
ment of optical thermal emission. For six of the eleven planets, 
there is little difference between the geometric albedo estimate af¬ 
ter a simple blackbody subtraction as opposed to the scenario with 
higher optical brightness temperature. For the remaining five plan¬ 
ets, the “full correction” geometric albedos are lower than their 
“simple correction” counterparts. The planets for which the details 
of thermal emission correction are more important tend to either 
have higher irradiation temperatures and so greater likelihood for 
unattributed thermal contamination in the optical (e.g. HAT-P-7b), 
or have precise optical eclipse measurements where minor changes 
to the dayside emission have a larger impact on constraining re¬ 
flected light (e.g. TrES-2b). Our contamination analysis is also 
largely consistent with the geometric albedos inferred from higher 
equilibrium temperatures in both Esteves et al.|j2014| and | Anger- 1 
|hausen et al.| ( |2014| >. Even in the unlikely event that hot Jupiters 
have optical dayside brightness temperatures equal to that in the 
mid-infrared, the optical geometric albedos for most planets are 
significantly lower than the Bond albedos inferred from thermal 
phase measurements. 

The second solution to the geometric versus Bond albedo dis¬ 
crepancy is that hot Jupiters are much brighter in the NIR, and 
hence thermal phase measurements—mostly obtained with Spitzer 
in the mid-IR—will underestimate their global temperatures and 
over-estimate their Bond albedos. However, neither the dayside 
emission spectra of individual planets, nor their aggregate spec¬ 
trum, show strong broadband molecular features. This means that 
dayside brightness temperatures from the near- through mid-IR 
should be reasonable proxies for their effective temperatures. One 
may worry that flux is escaping the nightsides of hot Jupiters in the 
NIR, leading us to underestimate nightside bolometric flux, but that 
is ruled out for WASP-43b by HST/WFC3 phase measurements, 
which show no nightside flux in the NIR (Stevenson et al.|2014| ). It 
would be useful to have full-orbit NIR phase curves of more planets 
in order to further test this hypothesis. 

This leaves us with the third hypothesis, namely that the Bond 


albedos of most hot Jupiters are high, despite their low geometric 
albedos. The geometric albedo of a planet (light reflected back to¬ 
wards the illuminating star) is related to its spherical albedo (light 
reflected in all directions) by a phase integral, A s = qA g . Lam¬ 
bertian (diffuse) reflection results in q = while pure Rayleigh 
scattering produces q = |. In general, planets with atmospheres— 
including simulated hot Jupiters—have 1.0 < q < 1.5 (jPollack] 
|et al.| 1986} [Burrows & Orton|201Q| l. It would be useful to empir¬ 
ically constrain the scattering phase functions of hot Jupiters us¬ 
ing data from space-based photometric missions. In the few cases 
where reflected phase variations have been measured, the spherical 
albedo appears so inhomogeneous that it is impossible to infer the 
phase-dependence of scattering (Demory et al. 2013, Esteves et al. 
|2014[ >. If we assume that hot Jupiters are diffusely reflecting, then 
they have typical optical spherical albedos of 15 per cent, still well 
below the inferred Bond albedos. 

The spherical albedo is related to the Bond albedo via a flux- 
weighted integral (jBurrows & Qrton|2010|l: 


/ 0 °°A s (A)/ lnc dA 
S JT^ncdA 


(5) 


where I mc is the SED of the incident stellar flux. The degree to 
which the optical spherical albedo impacts the Bond albedo de¬ 
pends on /“ , the percentage of starlight emitted in the observed 
optical waveband, assuming blackbody radiation at T* (values of 
/° pt are listed in Table[4]). We consider limiting cases of Equation^ 
assuming out-of-band wavelengths have average spherical albedos 
equal to 0, the optical A a , or 0.5 respectively: 


Ab 


A s f° 


opt 


(6) 


A gray _ a 

s 


(V) 


A“ 8h = A s fT + 0.5(1 - fT). 


( 8 ) 


Our limiting Bond albedos are summarized in Table[4] In principle 
the out-of-band spherical albedo could be unity, but as this would 
result in Bond albedos so great that the planets would be cooler than 
is observed in the mid-IR, we adopt a more modest upper threshold 
in Equation |8j 

The scenario is a reasonable match to the Bond albe¬ 
dos inferred from full-orbit thermal measurements. This suggests 
that most hot Jupiters have geometric albedos of « 50 per cent in 
the NIR and <10 per cent in the optical. If hot Jupiters are Lam¬ 
bertian reflectors, the NIR/optical contrast is somewhat less severe. 
This scenario similar in spirit to the high Bond albedo combined 
with low red-NIR geometric albedo one can obtain with Rayleigh 
scattering i jMarley et al.|1999| , but with the opposite color. Note 
that the high infrared albedos we are hypothesizing contradict the 
low infrared geometric al bedo we assumed when estimating re¬ 
flected IR light in Section . 


3.1.2 


Using A g = 0.5 in our MC im¬ 
plies greater NIR contamination from reflected starlight, and hence 
lower dayside thermal flux with greater Bond albedo. The most ex¬ 
treme change is a 20 per cent increase in the Bond albedo of WASP- 
12 b, but nonetheless our conclusions remain unaffected. 

It is worth mentioning that geometric albedos of 60 per cent— 
from the optical through the NIR—were predicted for the hottest 
giant exoplanets due to reflective silicate clouds l (Sudarsky et al.| 
|2000j ). In order to explain the low optical geometric albedo, one 
could invoke an optical absorber at low pressures, above the pur¬ 
ported cloud deck. Such optical absorbers, originally theorized to 
explain hot Jupiter temperature inversions, could include gaseous 


© 2015 RAS, MNRAS 000.|T|[T2] 




































Balancing the Energy Budget 11 


TiO/VO jFortney et al.|2007] l or S 2 /HS ( |Zahnle et al.|2009) . In this 
scenario, Kepler-7b would be unique not because of its clouds, but 
due to its dearth of optical absorbers. 

Alternatively, since cloud reflection is a multiple-scattering 
process, single-scattering albedos even marginally below unity can 
result in a heavily muted geometric albedo i jDlugach & Yanovitskij| 
|1974}|Hu|2014[ >. One might therefore explain the unusually red re¬ 
flectance spectrum of hot Jupiters with a single cloud deck where 
individual cloud grains have nearly gray albedo. 

Such scenarios might be tested by measuring the optical- 
infrared transit spectra of hot Jupiters. If the purported absorbers 
are located at sufficiently low pressures, and if the upper atmo¬ 
spheres of these planets are not too hazy ( |Pont et al. 2013, Gibson| 
|et al.|2013) , then we would expect larger effective radii in the opti¬ 
cal than the infrared. If instead the red reflectance spectrum is due 
to multiple scattering within a single cloud deck, then the transit 
spectrum should be flat. 

Moreover, the best way to investigate these hypotheses would 
be to obtain thermal phase measurements for the planets that have 
precise optical geometric albedo constraints, and vice versa. 
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